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Abstract 



Einstein's field equations for stationary Bianchi type II models with a perfect fluid source are inves- 



tigated. The field equations are rewritten as a system of autonomous first order differential equations. 
Dimensionless variables are subsequently introduced for which the reduced phase space is compact. The 
system is then studied qualitatively using the theory of dynamical systems. It is shown that the locally 
rotationally symmetric models are not asymptotically self-similar for small values of the independent 
variable. A new exact solution is also given. 
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I Introduction 

Dynamical systems methods have been used for over 20 years for studying the behavior of different models 
in general relativity, especially in the field of cosmology, see e.g. [El. The dynamical systems approach 
constitute a powerful tool when one wants to study asymptotic and intermediate behavior. It allows one 
to obtain a good understanding of the models even though it may be impossible to solve the corresponding 
equations exactly. Most of the attractors in the cosmological context have turned out to be self-similar 
solutions but there are also more exotic ones, e.g. the Mixmaster attractor. In this article we will apply 
the dynamical systems approach to the stationary Bianchi type II models. It will be shown that the locally 
rotationally symmetric (LRS) models are not asymptotically self-similar for small values of the independent 
variable. Instead the attractor is described by a heteroclinic cycle. To our knowledge, the stationary LRS 
type II models yield the simplest example of a non-self-similar attractor in general relativity. 

We will consider Bianchi type II models which admit a simply transitive group of isometries acting on 
3-dimensional hypersurfaces which are timelike. The line clement can be written as ds 2 — rj a bO} a U} b where 
the w a 's are the 1-forms 

o>° = Di{x) (dt + cydz) , U) 1 = dx , u 2 = D 2 (x)dy , u 3 = D 3 {x)dz , (1) 

and rjaf, = diag(— 1, 1, 1, 1). The parameter c is a constant. The source is assumed to be a perfect fluid for 
which the energy-momentum tensor has the form T a (, = ^u a u\, +p(i] a i, + u a Ub) where \i is the energy-density, 
p the pressure, and u a the 4- velocity of the fluid. The components of the fluid 4- velocity are u a = (1, 0, 0, 0). 
An equation of state of the form p(u) = (7 — l)/i with 1 < 7 < 2 is also assumed. The metric coefficients 
Dx, D 2 and D 3 are closely related to the kinematical properties of the normal congruence of the symmetry 
surfaces. Note that this congruence is spacelike. The expansion 9 and the shear a± are given by 

9 = -^ (In DxD 2 D 3 ) , 

dx 

+ 2 dx I D 2 D 3 J ' 2 dx I D 3 J K ' 



We also define 

Einstein's equations, G a b — T a b, lead to 
Evolution equations 
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Defining equation for /j 
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-|(6» + 4cr + )n , 
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-1)^= \9 2 - \o\ - \o 2 _ - n . 
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These equations can also be found from the orthonormal frame approach by specializing the equations of 
0] to the present models. The above set of equations is invariant under the transformation (9, a + ,a-,n) — > 
(9, cr+, — er_ , n). Therefore, without loss of generality, one can assume cr_ > 0. We can use the fact that fi is 
non-negative together with equation (pi) to see that 9 is a "dominant" quantity. Note also that 9, because 
of equation (||), cannot change sign. 

We now introduce 0-normalized variables: 

*-?. »-$. «-|f- W 



The introduction of a dimensionless independent variable r\ according to Odx = Zdr], leads to a decoupling 
of the ^-equation, 

9' = — = -(l + q)0 , q:=2Y? + + 2Y?_ + \{2- 1 )n . (10) 

The remaining equations can now be written in dimensionless form: 



Evolution equations 



Denning equation for Q, 



V + = -(2-g)£++4iV-7ft , (11) 

E_ = -(2-«)E_, (12) 

N' = 2{q-AT l+ )N . (13) 

(7- 1)0 = 1 -Si -S?. -N. (14) 



The boundary consists of a number of invariant sets which are important in understanding the dynamics 
of interior orbits and we will therefore include them. Moreover, this yields a compact reduced phase space. 
The boundary is given by (i) the static Bianchi type I models, N = 0, and (ii) the vacuum submanifold, 
fi = 0. There is also the locally rotationally symmetric (LRS) submanifold given by £_ = which divides 
the phase space into two parts, related by the discrete symmetry £_ — ► — £_. Note that the rotation of the 
fluid is non-zero, since lu 2 = n = N6 2 /3 H. 

II Dynamical systems analysis 

We start by listing the equilibrium points and the corresponding value of 0, which shows if the point is 
located on the vacuum submanifold or not. The eigenvalues for each point are also given but we refrain from 
giving the eigenvectors explicitly. 



The equilibrium points K 



Y, 2 + + Y, 2 _ = l , N = 0:n = (15) 

(57 ' 6)+ 1 27S+ , 4(1-22+), 0. (16) 
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The equilibrium point W 



3(5 7 2 ~ 4) _ = 3(7 7 -6) 
2(5 7 - 4) ' ' 4(5 7 - 4) 2 ' " (5 7 - A) 2 



E + = ^-^>E_=0,iV = ^^-^;O = ^^^ (17) 



12(7- I) 2 ±3v/A 12 ^-1 (18) 



2(57-4)(7-l) 57-4 

where 

A = 2( 7 -l)[8(7-l) 3 -(77-6)(5 7 2 -4)] . (19) 

The first two eigenvalues of W are always complex with a negative real part for 1 < 7 < 2. The point lies 
in the LRS submanifold and corresponds to the self-similar solution in |3| . 

Secondly we note that the equations corresponding to the static Bianchi type I boundary, N = 0, can be 
solved exactly. We find that 



where s± are constants satisfying s+ 2 + s_ 2 = 1. The vacuum boundary, tt = 0, is also solvable: 

£_ = S ~_ (£+ - 2) , N = 1 - Y? + - E 2 . . (21) 

A non-LRS exact solution which is characterized by a constant value of S + = (2 — 7)/(2(57 — 4)) can also 
be found (the solution was found using the Hamiltonian approach developed in [EJ. See the appendix for 
the explicit line element). The orbit starts on the Kasner circle K and ends at the point W . This solution 
is important as the orbits in the interior non-LRS part of the phase space spiral around it. The situation is 
analogous to that of the spatially homogeneous Bianchi type II models [g . 

The system of differential equations admits an increasing monotone function Z in the interior phase 
space, excluding the point W where Z takes its maximum value (this monotone function has been found by 
using the Hamiltonian methods developed in chapter 10 in [II]). The function is given by 

N m &- m 2-7 5 7 2 - 4 

(1- V S + ) 2 ' 2(5 7 -4)' (ll 7 -10)(37-2) ' l ' 

with 

^ = 1fi , (2(57-4)E + -(2-7)) 2 + 3(3 7 -2)(ll7-10)£ 2 _ 

Z {1 > (ll 7 -10)(3 7 -2)[2(57-4)-(2-7)E + ] " l > 

This monotone function prevents the existence of equilibrium points, periodic orbits, recurrent orbits and 
homoclinic orbits in this region, see e.g. J6|. 

The phaseportraits of the boundaries are given in figure la-d, while the phaseportrait of the LRS sub- 
manifold is shown in figure le. We see from the latter that orbits asymptotically approach W when rj — v oo, 
while for r/ — > — oo there exists a heteroclinic cycle, described by the LRS type I and the LRS vacuum 
submanifolds. Figure If depicts the phase space of the Bianchi type II non-LRS models. Note, in particular, 
the exact solution characterized by S + = (2 — j)/(2(5 r y — 4)). All other non-LRS orbits start at the Kasner 
circle, K, and spiral around this orbit towards W. In this case the LRS heteroclinic cycle is no longer an 
attractor. Instead it describes the intermediate behavior of those orbits which come "close" to it. 

So far we have only discussed the mathematical features of the stationary type II models. However, 
they might also be of some physical interest. Wainwright has speculated that the solution corresponding 
to the equilibrium point W might be interpreted as an approximation to the interior of a rotating disc of 
matter 0]. This interpretation should also pertain to the non-self-similar LRS models and perhaps also to 
the non-LRS models since they asymptotically approach W. However, the main importance of the present 
models is probably as part of a bigger picture where they may act as building blocks. The phase space of 
the present models forrnpart of the boundary of more general stationary Bianchi models and hypersurface 
self-similar models (see J2J). These models in turn form part of the boundary of more general models like the 
physically interesting G2-models, see e.g. 0. When studying these models one thus have to be observant of 
the behavior associated with the present heteroclinic cycle which could be expected to describe asymptotic 
or intermediate oscillating spatial behavior. 
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A The non-LRS exact solution 

The line element for the non-LRS exact solution is given by 

ds 2 = -[x(x + a)] 1/2 (dt + cydzf + [x{x + a)] q dx 2 

+x p + (x + a) p - dy 2 + x p - (x + a) p + dz 2 , (24) 
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Figure 1: The phase portraits of the stationary Bianchi type II models corresponding to (a) the N = 
boundary for 1 < 7 < 6/5, (b) the N — boundary for 7 = 6/5, (c) the N — boundary for 6/5 < 7 < 2, 
(d) the vacuum boundary, ft = projected onto the E±-plane, (e) the LRS submanifold, £_ — and (e) 
the full phase space. 



with 



37-4 
'4(7-1) 



P± 



V(ll7-10)(3 7 - 



37 - 2 ± S 
8(7-1) ' 
I „2 5 7 2 - 4 



16(7-1) 2 ' 



(25) 



where a is a constant. Setting a = yields the self-similar solution of |3|. For the remaining non-self-similar 
solutions one can set a = 1 by using the scale invariancc. 
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